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ABSTRACT 

We study a new s-process path through an isomer of 186 Re to improve a 
18 'Re- 18 'Os nucleo-cosmochronometer. The nucleus 18 'Re is produced by this 
new path of 185 Re(n, 7 ) 186 Re m (n, 7 ) 187 Re. We measure a ratio of neutron cap¬ 
ture cross-sections for the 185 Re(n, 7 ) 186 Re m and 185 Re(n, 7 ) 186 Re 9S reactions at 
thermal neutron energy because the ratio with the experimental uncertainty has 
not been reported. Using an activation method with reactor neutrons, we ob¬ 
tain the ratio of R t h = 0.54 ± 0.11%. From this ratio we estimate the ratio of 
Maxwellian averaged cross sections in a typical s-process environment at kT = 
30 keV with a help of the temperature dependence given in a statistical-model 
calculation because the energy dependence of the isomer/ground ratio is smaller 
than the absolute neutron capture cross-section. The ratio at kT =30 keV is esti¬ 
mated to be R st — 1.3 ± 0.8%. We calculate the s-process contribution from the 
new path in a steady-flow model. The additional abundance of 18 'Re through 
this path is estimated to be N s = 0.56 ± 0.35% relative to the abundance of 
186 0s. This additional increase of 18 'Re does not make any remarkable change in 
the 18 'Re- 18 'Os chronometer for an age estimate of a primitive meteorite, which 
has recently been found to be affected strongly by a single supernova r-process 
episode. 
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1. Introduction 

Two neutron-capture processes are important for astrophysical nucleosynthesis of heavy 
elements. The first one is a rapid neutron-capture process (r-process) that is considered to 
occur in supernova (SN) explosions (Meyer et ah 1992; Woosely et ah 1994; Qian, Vo¬ 
gel & Wasserburg 1998; Cowan et ah 1999; Otsuki, Tagoshi & Kajino 2000), and the 
other is a slow neutron-capture process (s-process) in low-mass asymptotic giant branch 
(AGB) stars (Kappeler et ah 1990; Straniero et ah 1995; Gallino et ah 1998; Arlandini et 
ah 1999) or massive stars (Prantzos et ah 1990; Woosley & Weaver 1995; Lee, ElEid & 
Meyer 2000; Hoffman, Woosely & Weaver 2001). Long-lived radioactive nuclei are used as 
nucleo-cosmochronometers, which are useful for an investigation of nucleosynthesis process 
history along the Galactic chemical evolution (GCE) before the solar system formation. A 
general idea of the nucleo-cosmochronometer was proposed by Rutherford about 70 years 
ago (Rutherford 1929). Radioactive nuclei of cosmological significance are rare and only six 
chronometers with half-lives in the range of the cosmic age 1 ~ 100 Gyr are known. They 
are 40 K (BBFH 1957) and 87 Rb (Clayton 1964) for the s-process or explosive nucleosynthesis 
in SNe, 176 Lu (Hayakawa et ah 2004) for the photodisintegration reaction nucleosynthesis in 
SNe (y-process), and 187 Re (Clayton 1964), 232 Th and 238 U (BBFH 1957; Fowler & Hoyle 
1960) for the r-process. Recently two actinoid elements LI and Th in a very metal-poor 
star were detected for the first time (Cayrel et ah 2001). The actinoid nuclei in the old 
metal-poor star were perhaps created in a single r-process of a SN explosion. Since the 
initial abundances of 238 U and 232 Th cannot be measured directly these abundances should 
be calculated in the framework of the r-process models. However the calculated results by 
many models (Otsuki, Tagoshi & Kajino 2000; Goriely & Arnould 2001; Schaty et ah 2002; 
Wanajo et ah 2002) are different from one another and thus a large uncertainty still remains 
in the estimation with the LI-Th chronometer. 

A pair of nuclei, 187 Re- 187 0s, was proposed to be an intentional cosmo chronometer to 
date the r-process (Clayton 1964). Fig. 1 shows a nuclear chart and flows of nucleosynthesis 
processes around the nucleus 187 Re. The nucleus 187 Re is predominantly produced by /3~- 
decay after the freezeout of the r-process and is located outside the main path of the s-process. 
The s-process nuclei 186 ' 18, Os are not directly produced by the r-process because they are 
shielded by the stable isobars of 186 W and 18 'Re against the /3-decay after the r-process. The 
pure s-process nucleus 186 0s is used for the normalization of the s-process abundances in this 



3 


mass region. The s-process nucleus 18 'Os is also produced by the cosmoradiogenic /3 _ -decay 
of 187 Re. The ground state of 18 'Re decays to 18 '0s with a half-life of 4.35xlO 10 yr (Lindner 
et al. 1986), which is longer than the age of the Universe. It should be noted that the 
18 'Re- 187 0s chronometer has the advantage that applying this chronometer is free from the 
uncertainty of the initial abundances calculated by the r-process models. The epoch from an 
r-process nucleosynthesis event to the present can be evaluated by the present abundances of 
18 'Re and 18 'Os after the subtraction of the s-process contributions to these nuclei. Although 
the initial abundances of 187 Re and 187 0s should be calculated using the s-process models, 
the uncertainty of the calculated s-process abundances is generally smaller than those in the 
r-process models. 

The 18 'Re- 187 0s chronometer has been applied for studying the GCE (Clayton 1969; 
Conser & Truran 1981; Yokoi, Takahashi & Arnould 1983; Arnould, Takahashi & Yokoi 
1984). This chronometer was also used in analyses of meteorites (Luck, Brich & Allegre 
1980; Luck & Allegre 1983; Birck & Allegre 1998). The measurements of the 187 Re and 187 0s 
abundances in iron meteorites and ordinary chondrites were used for evaluating the ages of 
the meteorites, which were formed from a metal reservoir at the solar system formation. 
These results lead to an upper-limit estimation of the age of the Galaxy (Luck, Brich & 
Allegre 1980; Luck & Allegre 1983). It was pointed out that the half-life of the /3 _ -decay of 
18 'Re is affected by stellar environments (Takahashi & Yokoi 1983). The half-life depends on 
the temperature and the electron density of the stellar environments. The half-life becomes 
shorter in high temperature such as T > 10 8 K due to the enhanced bound-state /3-decay in 
ionized atoms. This phenomena was verified in a measurement of the change of the half-life 
of the ionized 186 Re (Bosch et al. 1996). However the 18 'Re- 18 'Os system is still an important 
chronometer in the astronomical observations of old stars and in the analyses of pre-solar 
grains in primitive meteorites for the chronology of the r-process. 

The chemical compositions of the extremely metal-poor stars are enhanced by the pro¬ 
duction of heavy elements in neutron-capture nucleosynthesis in the first generation stars. 
Recently astronomical observations with isotope separations have been carried out for sev¬ 
eral heavy elements such as Ba (Lambert & Prieto 2002) and Eu (Sneden et al. 2002; Aoki 
et al. 2003a,b). Likewise, the pre-solar grains in the primitive meteorites provide samples 
affected strongly by a single nucleosynthesis event. Heavy elements such as Mo, Te and 
Ba in the pre-solar grains have already been separated into isotopes. The origins of these 
pre-solar grains are considered to be the ejecta of core collapse SN explosions (Richter, Ott 
& Begemann 1998; Pcllin et al. 2000; Yin, Jacobsen & Yamashita 2002). 

The purpose of this paper is to present a contribution of a new s-process path to 
18 'Re and 18 'Os in order to improve the accuracy of this cosmochronometer. In the pre- 
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vious studies the s-process contributions via two weak paths, 184 W(n, 7 ) 185 W(n, 7 ) 186 W(n, 7 ) 
187 W(z/,e _ ) 187 Re and 185 Re(n, 7 ) 186 Re ffS (n, 7 ) 187 Re, were studied. Kappeler et al. measured 
the neutron capture cross sections of 185 Re and 18 'Re at kT = 25 keV (Kappeler et al. 1991). 
The neutron capture cross section of an unstable nucleus 185 W is evaluated from the inverse 
reaction 186 W( 7 , n) 185 W (Sonnabend et al. 2003; Mohr et al. 2004). Here we propose another 
new weak s-process branch to produce 187 Re: a path through a 186 Re isomer which has a half- 
life of 2.0 xlO 5 yr. This isomer was found by Lindner about 50 years ago (Lindner 1951), 
and then Seegmillcr et al. measured the half-life of this isomer by a mass-spectrometric 
analysis and 7 -ray spectroscopy of samples after neutron irradiation (Seegmiller et al. 1972). 
However the absolute value of the neutron-capture cross section for the 185 Re(n, 7 ) 186 Re m 
reaction has not been reported. Seegmiller et al. reported a ratio of the neutron-capture 
cross section of 185 Re leading to the isomer and the ground state, but the uncertainty was 
not given in the article. Therefore the neutron capture cross section to the isomer has not 
been established (Firestone et al. 1998). The half-life of the ground state of 186 Re is only 
3.718 days, while that of the isomer, 2.0 xlO 5 yr, is longer than the timescale of the s-process, 
10 2 ~ 10 4 yr. The s-process contribution should be, thus, estimated with high accuracy. 

In the present work we measure the neutron-capture cross section ratio 185 Re(n, 7 ) 186 Re m 
/ 185 Re(n, 7 ) 186 Re 9S at the thermal energy. An implication of this path in the abundances of 
18 'Re and 18 'Os is discussed using our data in a classical steady-flow model, and an effect to 
the chronometry is evaluated in a sudden approximation. 


2. Experiment procedure and results 
2.1. Measurement of neutron capture cross-section ratio 

In order to measure the ratio of the neutron-capture cross-sections 185 Re(n, 7 ) 186 Re m 
and 185 Re(n, 7 ) 186 Re 9S at the thermal energy, two thin natural Re foils were irradiated with 
thermal neutrons from the nuclear reactor JRR-4 at the Japan Atomic Energy Research 
Institute (JAERI). Two Al-Co wires were also irradiated simultaneously to monitor the 
neutron flux. The samples were irradiated by the reactor neutrons for six hours. After the 
irradiation the samples were cooled in a water pool for about four months in order to reduce 
the background of the unstable isotopes with short half-lives. Four and eight months after 
the irradiation, 7 -rays emitted from the activated samples were measured by two HPGe 
detectors. In addition sixteen months later, the 7 -rays from one sample were measured by 
one HPGe detector. 

The weights of the Re foils are 43.8 mg and 58.1 mg, whereas the weights of the Al- 
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Co wires are 1.74 mg and 4.58 mg. The Co element composition is 0.475%. The natural 
Re sample consists of two stable isotopes, namely 185 Re (37.4%) and 18 'Re (62.6%), and 
thus two unstable isotopes, 186 Re and 188 Re, are mainly produced by the neutron capture 
reactions. The half-lives of the ground states of these isotopes are T { / 2 = 3.718 days ( 186 Re) 
and T \/2 = 17 hr ( 188 Re), which are much shorter than the 186 Re isomer. The average flux of 
the neutrons is estimated to be 4.37xl0 13 cm -2 s -1 , which is evaluated by a measurement of 
decay y-rays after /3-decay from the activated Co in the wires. The y-rays from 28 A1 cannot 
be measured because a cooling time of four months is much longer than the half-life of 28 Al. 
The total efficiency of the HPGe detectors is approximately equal to 20% relative to a 3” x3” 
Nal detector and the energy resolution is 2.2 keV for a 1.3 MeV y-ray in b0 Co. The efficiency 
is calibrated with standard sources of 152 Eu and 133 Ba. The y-rays from the samples are 
measured at a 5 cm or 3 cm distance from the end caps of the detectors in the first or second 
(third) measurement, respectively, with a lead-shielded low-background environment. 

Fig. 2 shows a partial decay scheme of the 186 Re isomer (Baglin 2003). The isomer 
with J 7r =(8 + ) locates at an excitation energy of 149 keV. This state decays to the ground 
state through a cascade of y-decay. The ground state subsequently decays to 186 0s or 186 W 
with 7\/2=3.718 days. The y-ray transitions from the isomer to the ground state cannot be 
distinguished from y-rays emitted from other activities. Strong y-rays are irradiated from 
184 Re (T 1 /2=38 days), which is produced by a (n,2n) reaction with fast neutrons existing in 
the nuclear reactor (see Fig. 3). Its half-life is comparable to the cooling time, and hence 
the intensities are enhanced. However we observe a 137 keV transition of 186 0s (see Fig. 4), 
which is the strongest y-ray in the decay scheme associated with the isomer. The yield of 
the 137-keV y-ray is contributed not only from the isomer but also from the ground state 
produced directly by the l85 R,e(n,y) 186 Re 9S reaction. 

In order to obtain the ratio of the neutron-capture cross section between the isomer 
and the ground state, we select the cooling time of four months after the irradiation for the 
first measurement. After the cooling, the ground state of 186 Re is still remained because of 
its large capture cross-section relative to the isomer. Although the surviving number of the 
ground state is small, the contribution to the 137 keV y-ray yield is competitive to the isomer 
because the decay rate of the ground state is much higher than the isomer. In the second and 
third measurements the ground states are not remained. The yields of the 137-keV y-ray 
at the second and third measurements are, thus, proportional to the neutron capture cross 
section to the isomer. 

The number of the nuclei produced by the neutron-capture reaction is proportional to 
the /3-decay events, which are evaluated from the y-ray intensity in the measured spectra by 
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applying the following equation of 


Ntotal 


Tlpeak 


R-e-f dl 


ecay 


(i) 


where N tota i means the total number of the /5-decays, n pea k denotes the peak counts corre¬ 
sponding to the decay 7 -rays, R the emission probability of the measured 7 -ray per decay, 
e the efficiency of the HPGe detectors and fdecay the correction of the /5-decay during the 
measurement time. The fdecay is expressed by 


fdecay = exp(-5 r A) - exp(-t 2 -A), 


( 2 ) 


where t\ and t 2 stand for the start and stop time of the measurement, respectively. The 
experimental uncertainty consists mainly of the statistical error of the yields of the 7 -rays 
and the efficiency of the HPGe detectors. The 7 -ray energy of 137 keV is near the peak of the 
efficiency curve of the detectors. We evaluate the efficiency uncertainty to be approximately 
7% from the deviation of the 7 -ray efficiency obtained with the calibration sources from the 
efficiency curve estimated by y 2 fitting. The pileup effect of the 7 -rays is negligibly small 
because the absolute efficiency is lower than 5%. The self-absorption of the 7 -rays is also 
negligibly small. 


An energy spectrum of neutrons produced by a nuclear reactor consists of thermal flux 
with energy below 0.2 eV and epithermal one above. The rate of the neutron capture reaction 
is expressed by 

N x = J n(v)va(v) x dv = (j) th crQ + , (3) 

where N x denotes the reaction rate for the reaction x, n{y) the neutron density at the velocity 
v, <j(v) x the cross section of the reaction x at the velocity v, dth and 4> ep i the thermal and 
epithermal flux, aft and Jq are the cross section at the thermal energy and the resonance 
integral for the reaction x, respectively. A measurement for a sample with a Cd absorber is 
not carried out in the present experiment and we cannot obtain the thermal neutron flux. 
Instead we use an average neutron flux, which is defined by 


_ Ntotai ( 60 Co) 
W “ M 60 Co) 


( 4 ) 


where (j) ave the average neutron flux and a x ^ the effective cross section of the reaction x in 
the reactor neutron spectrum. The effective neutron-capture cross sections with the average 
neutron flux are defined by 
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Applying these equations, we obtain the effective neutron capture cross sections leading to 
the isomer of 186 Re as given in Table 1, where the mean value derived from the second and 
third measurements is 0.74 ± 0.05 barn. The yield of the y-ray at the first measurement 
consists of contributions from the ground-state decay as well as the isomer decay. The 
effective neutron capture cross section to the ground state can be evaluated by subtracting 
the contribution of the isomer decay, which is calculated from the average cross section of 
0.74 ± 0.05 barn, from the total yield at the first measurement (see Table 2). The obtained 
cross section to the ground state is 132 ± 26 barn. This value is consistent with the thermal 
values of 118 barn (Nakajima 1991) and 112 barn (Mughabghab 1984) within the uncertainty 
and thus the effective cross section derived here is essentially the thermal cross section. 
To verify this conjecture, we compare the thermal cross sections and resonance integrals 
leading to the ground state and isomer of 186 Re and their ratio in Table 3. These values 
are calculated from Japanese Evaluated Nuclear Data Library (JENDL) Activation Cross 
Section File 96 (Nakajima 1991; JENDL 2000). The ratios, Jo/oo, for both the ground state 
and isomer are approximately equal (see Table 3). This fact indicates that the resonance 
integral are written by 


= aa o S , 

I™ = aa™, 


(7) 

( 8 ) 


with the common coefficient a. The isomer/ground ratio of the capture cross sections at the 
thermal energy can be derived as 

g/T _ <7o l (&th + a<f>epi) _ V™<Pth + 7p % <t>epj _ , q x 

cr 0 * Vo S {<t>th + OKpepi) CTo* <j>th + Io S <f>epi of//’ 

which verifies our conjecture. The isomer/ground ratio obtained by applying this equation 
is R t h = 0.54 ± 0.11%. This value is a factor of 1.8 larger than the value of 0.3% reported 

in a previous report (Seegmillcr et al. 1972). Note that experimental uncertainty was not 

reported in this report (Seegmiller et al. 1972). 


There are strong y-rays from 184 Re produced by the 185 Re(n,2n) 184 Re reaction with 
fast neutrons in Fig. 3. The half-life of 184 Re (77/2=38 days) is long enough to remain 
after the cooling time. This fact suggests the possibility that the isomer of 186 Re might 
be also produced by the 187 Re(n,2n) 186m Re reaction. To estimate the contribution by the 
18 'Re(n,2n) 186 Re m reaction, we evaluate the effective cross section of the 185 Re(n,2n) 184 Re 


reaction by 

(n,2n) / 184p \ Re) 

a eff ( Re ) = —7 - 


( 10 ) 


We obtain the result that the effective cross section of the 185 Re(n,2n) 184 Re reaction is 
approximately equal to 0.00078 barn. Since the decay rate of 184 Re is much higher than 



186 Re m , the 7 -ray intensities of 184 Re are stronger than the 186 Re isomer. Assuming that 
the effective cross section of the 187 Re(n,2n) 186 Re is equal to that of the 185 Re(n,2n) 184 Re 
reaction, the contribution to the 186 Re isomer can be calculated by 

lV( 186 Re m ) = ( 11 ) 

where R means the production ratio of the isomer to the ground state of 186 Re, A(x) means 
the isotope abundance ratio in the sample and N ( x ) means the number of the nuclei produced 
by the neutron-induced reactions. The contribution of the 187 Re(n,2n) 186 Re m reaction is ap¬ 
proximately equal to 0.0007% relative to that of the 185 Re(n, 7 ) 186 Re m reaction. We therefore 
conclude that the contamination due to the 187 Re(n, 2 n) 186 Re TO reaction is negligibly small. 


3. Discussion 

3.1. s-Process branch to 18 'Re through the 186 Re isomer 

We here estimate the effect of the new path in the s-process using a classical steady-flow 
model (BBFH 1957; Clayton 1964). The steady-flow model was extended by Ward, Newman 
and Clayton (Ward, Newman & Clayton 1976) to apply to a branching point where the s- 
process nucleus is produced by a neutron capture and/or a feeding by /5-decay of a parent 
nucleus and also destroyed by a neutron capture and/or /5-decay. The abundance of the 
residual nucleus can be calculated by 

5^24 = Y,HA)N‘(A) - y>(V)A P(A'), ( 12 ) 

" * 3 

where A means the neutron capture rate ( n-a ■ v ) or the /5-decay rate ( ln2/Ti / 2 ), n is the 
neutron number density, T^ is the half life of the /5-decay, N means the isotope abundance, 
A and A' are the nuclear mass numbers which have the relationship, A=A' (for the /5-decay) 
or A'-l (for the neutron capture). In the limit that the timescale of the s-process is much 
longer than a mean time of neutron capture, dN/dt vanishes and hence the ” steady flow” 
becomes a good approximation of the s-process. The equation in this limit is expressed by 

5>e»)V(.4) = J2W)W(A'). (13) 

i 3 

Applying this equation to the Re-Os branching point displayed in Fig. 1, we obtain a set of 
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the following equations: 


A ( g-iV s ( 186 Re 9S ) = A(to 187 0s)-/V s ( 186 0s), (14) 

Xp-N s ( 1S5 W) = A (to 186 Re m ) AG ( 185 Re) + A (to 186 Re 9S )W( 185 Re), (15) 

A (to 186 Re 9S )Y s ( 185 Re) = (Xp- + A EC//3+ )lV s ( 186 Re 3S ), (16) 

A(to 186 Re m ) N s ( 185 Re) + A(to 187 W)W( 186 W) = A(to 188 Re)Y s ( 187 Re), (17) 

A(to 185 W)W( 184 W) = Xp-N s ( 1S5 W) + A(to 186 W)W( 185 W), (18) 

A / 3 +lV s ( 186 Re 9S ) + A(to 186 W)W( 185 W) = A (to 187 W)1V,( 186 W). (19) 


We here assume that the internal transition rate between the ground state and the isomer 
is negligibly lower than the neutron capture rate. The validity of this assumption should be 
discussed later. A typical astrophysical nucleosynthesis environment of the s-process is kT = 
30 keV and n = 10 8 cm -3 . Since the following two relations hold, A(to 188 Re) 3> A j a-( 18 'Re) 
at 187 Re and (A/j- + A ec//3+ )( 186 Re 3S ) 3> A(to 187 Re) at 186 Re 3S , on the environmental 
condition for T!/ 2 ( 187 Re) = 43.5 Gyr and T 1 / 2 ( 186 Re s ' s ) = 3.718 days, we omit those two 
slower reaction rates in Eqs. (16) and (17). The situation does not change when the atom 
of 18 'Re is ionized in the stellar interior during the s-process and the half-life decreases by 
9~10 orders of magnitude (Yokoi, Takahashi & Arnould 1983; Bosch et ah 1996). Two weak 
s-process contributions are known (Kappcler et al. 1991). The unstable nucleus 185 W is a 
branching point where it decays to 185 Re with a half-life of 75.1 days and it is transmuted to 
186 W by the neutron capture reaction. The nucleosynthesis flow from 186 W reaches to 18 'Re. 
The /3-decay rate at 185 W is higher than the neutron capture rate by an order in the case 
of n=10 8 cm -3 , but it depends on the neutron density in the s-process environment. The 
ground state of 186 Re is the second branching point, which decays to both 186 0s and 186 W. 
Using Eqs. (14)-(19), the s-process abundance of 187 Re relative to 186 0s that is produced by 
the s-process including the new s-process path is calculated by 


lY s ( 187 Re) 


R A g+TRen 1+R A(to 186 H-') A(fo 187 0s) 18e 

t J 4 A j3 -( ls6 /?eB>) A X(to lm Re) h 


( 20 ) 


where R is defined by A(to 186 Re m )/A(to 186 Re 9S ) and A is the branching ratio of the (3- 
decay of 186 Re 3 *, namely A = Xp-/(\p-+\ ec// 3 +)- We use the branching ratio of 0.9253 for 
A (Baglin 2003). It should be noted that the s-process abundance of 187 Re has two terms 
proportional to the ratio of the neutron capture rates leading to the isomer and the ground 
state of 186 Re. 
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The Maxwellian-average neutron-capture cross section (MACS) with quasi-stcllar neu¬ 
trons corresponding to an energy range of kT— 8~30 keV is important for calculating nu¬ 
cleosynthesis under typical s-process environments. We measure only the neutron-capture 
cross-sections with the thermal neutrons in the present experiment. In order to estimate 
the s-process abundance of 18 'Re, we assume first that the ratio of neutron capture rates, 
A(to 186 Re m )/A(to 186 Re 9S ) = 0.54 ± 0.11%, does not depend on the neutron energy. The 
MACS at kT —30 keV for Re and Os isotopes are taken from the best set evaluated pre¬ 
viously (Kappeler et al. 1991). The values are 1160 and 418 mbarn for 18 'Re and 186 0s, 
respectively. We assume the neutron density of 10 8 cm -3 to estimate a neutron capture rate 
of A(to 186 W). Substituting the measured cross section ratio, R t h = 0.54%, into Eq. 20, we 
obtain the s-process abundance of 18 'Re of N s = 0.23 ± 0.05% relative to 186 0s. 

As an alternative method to estimate the isomer/ground ratio at stellar energy, we 
calculate the energy dependence of the ratio by the Hauser-Feshbach statistical-model, which 
was used in JENDL Activation Cross Section File (Nakajima 1991; JENDL 2000). In JENDL 
the total capture cross section is calculated using a resonance formula and the statistical- 
model (see Fig. 5). We first calculate the capture cross section to the ground state, which is 
normalized at 100 keV to the value given in JENDL. The isomer/ground ratio is also given 
in JENDL, but its value at the thermal energy is smaller than the presently measured value 
by a factor of 23. In order to reproduce the measured isomer/ground ratio at the thermal 
energy, we adjust the parameters in the present statistical-model calculation in two different 
ways. In the first method we adjust the E 3 transition probability in the q-decay cascade 
to enhance the production of the isomer. In the second method we adjust the level density 
parameters (LDP) which affects the branching of the neutron and 7 emission and affects the 
average number of cascade 7 -transitions. The calculated isomer/ground ratios are shown in 
Fig. 6 as a function of the incident neutron energy. The fact that the energy dependence 
of the two calculated results is almost identical indicates that the calculated ratio is not 
sensitive to the choice of the parameters. Finally both results give the isomer/ground ratio 
of the MACS at kT =30 keV to be R st = 1.3 ± 0.8%. Although the two calculated ratios 
are identical, we assign an error of 50% to the calculated result as a conservative estimate of 
the uncertainty. By inserting this ratio, R st = 1.3 ± 0.8%, to Eq. 20 we obtain the s-process 
abundance of 18 'Re, N s = 0.56 ± 0.35%, relative to 186 0s. This value is larger than that at 
the thermal energy by a factor of 2.4. This difference shows that the energy dependence of 
the ratio is important for the study of the s-process and an experimental measurement of 
the cross section at the stellar energy is a further subject. 

Fig. 6 shows that the calculated ratio increases drastically above 10 keV. This tendency 
indicates that the abundance ratios for the Re and Os isotopes can be used as a nucleo- 
thermometer of the s-process due to the strong temperature dependence of the isomer/ground 
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ratio. 

Finally we would like to discuss a problem concerning the transition probability between 
the isomer and the ground state. We assume that this transition probability on a typical 
s-process condition is ignored. A case of the nucleus 180 Ta gives a clue for this problem. The 
nucleus 180 Ta has a similar nuclear structure: the isomer is meta-stable with J 7 r =9 + (7j / 2 > 
10 15 yr) existing in the solar system, while the ground state with J n =l + is unstable against 
the /3-decay with a half-life of 8.15 hr. This nucleus is proposed to be synthesized through 
a weak branch of the s-process (Yokoi & Takahashi 1983; Belie et al. 1999) and by photo¬ 
disintegration reactions in supernova explosions (Arnould, & Goriely 2003; Utsunomiya et 
al. 2003). The isomer synthesized by either process may decay to the ground state through 
intermediate states (IMS) located above the isomer by ( 7 , 7 ') reactions in high tempera¬ 
ture environment, and subsequently the ground state decays to daughter nuclei through 
/3-decay. The residual abundance of 180 Ta is affected by the transition probability between 
the ground state and the isomer, which depends on the excitation energy of the IMS (Belie 
et al. 1999). It is pointed out that the 1S0 Ta abundance contributed by the s-process should 
be smaller than the present solar abundance if the IMS might exist below excitation energy 
of 1 MeV (Lee, ElEid & Meyer 2000). This suggests that there is no IMS lower than 1 MeV. 
The IMS have been investigated experimentally, but the IMS below 1 MeV has not been 
found. Recently the transition probability from the isomer to the ground state has been 
directly measured in photo-excitation experiment (Belie et al. 1999) and the result supports 
the previous investigation. 

In the case of 186 Re, the IMS at low excitation energy also play an essential role for the 
destruction of the synthesized 186 Re isomer, but these IMS of 186 Re have not been found. 
The transition probability between the ground state and isomer via IMS depends on the 
quantum number K (Bear 1981). The quantum number K is a projected component of 
the total angular momentum on the nuclear symmetry axis. The transition probability 
increases with increasing the deformation of the nuclear shape (Carroll et al. 1991), because 
the Coriolis interaction in high spin states of the deformed nuclei causes the mixing of the 
states with the different K values. Since 180 Ta is more distant from the magic number, 
the deformation of 180 Ta is larger than 186 Re (Nazarewicz 1990). This fact indicates that 
the effect of the Coriolis interaction on 186 Re is smaller than that of 180 Ta. We therefore 
assume that the transition between the ground state and the isomer via the IMS can be 
ignored in the typical s-process environment. The measurement of the transition probability 
of 186 Re may be important for the study of the s-process abundance of 18 'Re. The direct 
measurement of the transition probability is difficult because the isomer of 186 Re does not 
exist in the solar system. A search for excited states by in-beam 7 -ray spectroscopy such as 
a measurement of excited states of 18 'Re (Shizuma et al. 2003) would be a further subject. 
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3.2. Effect on 187 Re- 187 Os chronometry 

The 18 'Re - 187 0s chronometer is useful for the study of the age estimate of an r-process 
event. The purpose of this subsection is to discuss the effect of the new s-process path on the 
18 'Re - 187 0s chronology. The contribution from the new s-process path through the 186 Re 
isomer has been totally ignored in the past. We first present the age estimate without the 
effects of all weak s-process branchings. In such a case an age of an object which exhibits 
both the 187 Re and 187 0s abundances is calculated by 

r 1 / 2 ( 187 Re) 

In 2 

AU 187 Re ) 

AU 187 0s) - N s ( 187 0s) + AU 187 Re) 

where N ob (A ) means the observed abundance and Ns (A) means the abundance contributed 
from the s-process. Note that we here adopt a simple model taking a sudden approximation 
that the r-process occurred at the look back time at t — T only once. Although the sudden 
approximation is not generally a good approximation for analyses of the solar materials, 
it is still useful in analyses of primitive meteorites that have recently been found to be 
affected strongly by a single supernova r-process episode (Amari et al. 1992; Yin, Jacobsen 
& Yamashita 2002). This approximation also applies to an analysis of astronomical date 
such as the abundances of the r-process elements detected in metal-poor stars with isotope 
separation because the elements on the surface of the old low-metallicity stars are probably 
generated in a single nucleosynthesis event (Sneden et al. 1996; Cayrel et al. 2001; Otsuki, 
Mathews & Kajino 2003; Honda et al. 2004). 

We present the age estimation by taking into account the new s-process path to 18 'Re. 
In this case, Eq. (21) should be replaced by 




T n/2( 187 Rey 

ln2 

( ^( 187 R«)-^s( 187 Re)-exp(-T s -ln2/T 1/2 ( 187 ite)) \ 
n VAU( 187 0s) - iV s ( 187 0s) + (AU( 187 Re) - Y s ( 187 Re)) )’ 1 j 

where we assume that the s-process occurred at different look back time T$ from the r- 
process episode, i.e. Tj^T s , N S (A) means the contribution by the s-process, N ob (A ) means 
the observed abundance. The s-process contribution to 18 'Os should be modified by the 
change of the nucleosynthesis flow at 186 Re. The s-process contribution to 18 'Re presented 
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in Eq. (20) is calculated from Eqs. (14) ~ (19), which contain all the weak s-process paths. 
In Eq. (20) the contribution of the new s-process path is expressed by the term of R. 

Substituting the estimated MACS ratio of R st = 1.3% at kT= 30 keV, we obtain the 
results shown in Tables 4 and 5. The neutron capture cross sections at kT = 30 keV 
are taken from the previous article (Kappcler et al. 1991). Those of 186 ’ 18, Os and 18 'Re 
are 418, 874, 1160 mbarn, respectively. We use the solar abundances for Wfe( 18 'Re) and 
N ob ( 187 0s). A(s( 187 0s) is calculated from the observed abundance Wfe( 186 Os) = A(s( 186 0s). 
The Wb( 186 0s)/lV o f>( 18 'Re) ratios in the primitive meteorites or the metal-poor stars may be 
different from the solar abundance ratio because 7V 0 &( 187 Re) should be enhanced in the sample 
which is strongly affected by the single r-process event. We thus introduce a parameter, /, 
which is defined by 


A r o b( 186 0s)/A r o fe( 187 Re) = /xlV 0 ( 186 Os)/lV 0 ( 187 Re), (23) 

for allowing the variation of N ob ( 186 0s)/N ob ( 187 Re). 

The difference between the calculated ages for R = 0 and 1.3% is at most ~ 1% in 
either case Ts=5xl0 9 yr (Table 4) or 10 10 yr (Table 5). The ages are slightly sensitive to 
the neutron density. Aoki et al. reported isotope ratios of Eu in s-process element rich 
metal-poor stars and show the possible neutron density range of the s-process environment 
by comparing their data with the s-process calculation: 10' < n < 10 9 cm -3 (Aoki et al. 
2003a,b). Thus we also present the age for the neutron density of 10 7 cnr 3 as well as 10 s 
cm -3 in Tables 4 and 5. The difference between the ages for two different neutron densities 
is at most ~ 2%. Several metal-poor halo stars are thought to form from either the s-process 
element rich gas (Aoki et al. 2003a) or the r-process element rich material (Sneden et al. 
1996; Cayrel et al. 2001; Honda et al. 2004), which are ejected from progenitor stars. The 
parameter / can be even larger or smaller than unity. The primitive meteorites, which are 
affected by a single supernova r-process, can have the / values smaller than unity. The 
column of / in Tables 4 and 5 shows this parameter. The estimated ages for R — 0 and 
1.3% under the conditions of / < 1.0 show almost the same results because the absolute 
nucleosynthesis flow from the s-process decreases with decreasing the / value. In particular, 
the ages become slightly different with increasing / value for the case of Ts=10 10 yr (Table 
5). Although the effect of the new s-process path seems to be small, the effect is almost 
the same as those of the other s-process paths under the condition of / < 1.0. Therefore 
the direct measurement of the absolute neutron capture cross section of 185 Re leading to the 
isomer of 186 Re at kT= 8~30 keV is important for the 18 'Re - 18 'Os chronometer. 



14 


4. Summary 

We propose a new s-process path that leads to synthesize 18 'Re and 18 'Os, which form 
an important nuclear-pair as a nucleo-cosmochronometer of the r-process. The new path 
consists of a neutron capture reaction chain from 185 Re to the isomer of 186 Re and from the 
186 Re isomer to 18 'Re, namely 185 Re(n, 7 ) 186 Re m (n, 7 ) 18 'Re, which has been ignored in all 
previous studies. We measure the ratio of the reaction cross sections 185 Re(n, 7 ) 186 Re m and 
185 Re(n, 7 ) 186 Re 9S with thermal neutrons provided by a nuclear reactor at the JAERI. We 
measure 7 -rays after the /3-decay from the Re samples by HPGe detectors and thereby we 
obtain the ratio of Rth = 0.54 ± 0.11%. The Maxwellian averaged cross section reaction 
ratio at the stellar energy kT = 30 keV is estimated with the help of the statistical model 
calculation. The result is used to calculate the s-process abundance of 187 Re relative to the 
abundance of 186 0s in the classical steady-state flow model, which yielded a value of N s = 
0.56 ± 0.35%. The proposed new s-process path does not make any remarkable change in 
the 18 'Re- 18 'Os chronometer in the sudden approximation when all other possible s-process 
branchings are included in the calculation. This confirms the robustness of the 187 Re- 187 0s 
chronometer which applies to the age estimated by the analyses of primitive meteorites that 
are presumed to be affected strongly by a single supernova r-process episode. Since the 18 'Re 
atom in these meteorites most likely has not been astrated in the stellar interior after it was 
produced in the r-process, the chronometer is free from the change of the half-life of ionized 
18 'Re. We however need further experimental study to measure the cross section precisely 
at the stellar energy kT— 8~30 keV. 

We would like to thank K. Takahashi, H. Utsunomiya and M. Fujiwara for valuable 
discussions. We also thank the crew of the nuclear reactor at JAERI. This work has been 
supported in part by Grants-in-Aid for Scientific Research (15740168) of the Ministry of Ed¬ 
ucation, Culture, Sports, Science and Technology of Japan, and the Mitsubishi Foundation. 
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Table 1: Effective neutron-capture cross sections of 185 Re leading to the 186 Re isomer. Two 
samples were irradiated for about six hours by reactor neutrons. Time means the cooling 
time after the irradiation. Since the cooling time is long enough, the ground state with a 
short half-life does not survive._ 


Sample 

Time 

(month) 

(barn) 

error (barn) 

1 

8 

0.78 

0.06 


12 

0.60 

0.08 

2 

8 

0.84 

0.10 

Average 


0.74 

0.05 
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Table 2: Effective neutron capture cross section of 185 Re leading to the ground state of 
186 Re. Two samples were irradiated for about 6 hours by reactor neutrons. Time means the 
cooling time after the irradiation. The capture cross section was obtained by subtracting 
the contribution of the isomer decay. 


Sample 

Time 

% 7 / (barn) 

error (barn) 


(month) 



1 

4 

136 

30 

2 

4 

142 

52 

Average 


132 

26 


Table 3: Neutron capture cross-sections at thermal energy, the resonance integrals and their 
ratios for the 185 Re(n, 7 ) 186 Re 9S and 185 Re(n, 7 ) 186 Re m reactions calculated from JENDL Ac- 


tivation Cross Section File 96 


a 0 (barn) 

J 0 (barn) 


gs 117.6 

1850 

15.7 

isomer 0.02711 

0.443 

16.3 


Table 4: Comparison of the ages in units of Gyr calculated using the 18 'Re - 18 'Os chronometer 
with (R = 1.3%) and without (R = 0%) the new s-process path. The age of the s-process, 
T s = 5 xlO 9 yr, is assumed. / is a parameter to allow the variation of observed 186 0s/ 18 'Re 
ratio relative to the solar ratio (see Eq. (23) ) and n is the neutron density in units of cm -3 . 


R means the ratio of the neutron capture leading the isomer and the ground state of 186 Re. 


/ 

Eq. (21) 

Eq. (22), n=10 8 , R=0% 

Eq. (22), n=10 8 , R=1.3% 

Eq. (22), n=10 7 , R=1.3% 

1.5 

4.86 

4.86 

4.86 

4.86 

1.2 

7.73 

7.80 

7.80 

7.76 

1.0 

9.56 

9.66 

9.67 

9.61 

0.5 

13.9 

14.0 

14.1 

14.0 

0.1 

17.2 

17.3 

17.3 

17.2 


Table 5: The same as those in Table 4, but for Ts=10 10 yr 


/ 

Eq. (21) 

Eq. (22), n=10 8 , R=0% 

Eq. (22), n= 10 8 , R=1.3% 

Eq. (22), n=10 7 , R=1.3% 

1.5 

4.86 

4.69 

4.68 

4.79 

1.2 

7.73 

7.67 

7.66 

7.70 

1.0 

9.56 

9.55 

9.55 

9.56 

0.5 

13.9 

14.0 

14.0 

14.0 

0.1 

17.2 

17.2 

17.2 

17.2 
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\ \ \ \ \ 


r-process 


Fig. 1.— Nuclear chart and nucleosynthesis flow around 18 'Re. The nuclei in this mass 
region are synthesized by the s- and r-processes. The ground state of 18 'Re decays to 18 'Os 
with a half-life of about 44 Gyr. The dotted-line is the new s-process path to 18 'Re proposed 
in this work. 
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Fig. 2.— Partial decay scheme of the 186 Re isomer. Relative 7 -ray intensities are indicated. 
Large number means the energy of the 7 -ray and small number means the relative inten¬ 
sity. The isomer of 186 Re decays to the ground state through the internal transitions and 
subsequently the ground state disintegrates to 186 0s or 186 W. 



Counts per Channel g; &■ Counts per Channel 


23 



0 200 400 600 800 


Energy [keV] 


3.— Spectrum of 7 -rays from the nat Re target. Strong 7 -rays from 184 Re are observed, 
h is produced by the 185 Re(n, 2 n) 184 Re reaction. 



Fig. 4.— The 7 -ray of 186 0s. The 137 keV 7 -ray is clearly observed. 
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Fig. 5.— Evaluated neutron capture cross section data for the 185 Re(n, 7 ) reaction leading 
to the ground state (solid line) and isomer (dashed line) of 186 Re given in JENDL Activation 
Cross Section File 96. 
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Fig. 6.— Cross section ratios of the isomer to the ground state of 186 Re calculated by the 
Hauser-Feshbach statistical-model. These ratios are adjusted to reproduce the presently 
measured value at the thermal neutron energy (see text). 




